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The evolution and population dynamics of avian coronaviruses (AvCoVs) remain underexplored. In the 
present study, in-depth phylogenetic and Bayesian phylogeographic studies were conducted to investi¬ 
gate the evolutionary dynamics of AvCoVs detected in wild and synanthropic birds. A total of 500 samples, 
including tracheal and cloacal swabs collected from 312 wild birds belonging to 42 species, were ana¬ 
lysed using molecular assays. A total of 65 samples (13%) from 22 bird species were positive for AvCoV. 
Molecular evolution analyses revealed that the sequences from samples collected in Brazil did not cluster 
with any of the AvCoV SI gene sequences deposited in the GenBank database. Bayesian framework anal¬ 
ysis estimated an AvCoV strain from Sweden (1999) as the most recent common ancestor of the AvCoVs 
detected in this study. Furthermore, the analysis inferred an increase in the AvCoV dynamic demographic 
population in different wild and synanthropic bird species, suggesting that birds may be potential new 
hosts responsible for spreading this virus. 

© 2015 Elsevier B.V. All rights reserved. 


1. Introduction 

Avian coronaviruses (AvCoVs) are enveloped viruses with 
pleomorphic spherical forms belonging to the order Nidovirales, 
family Coronaviridae, subfamily Coronavirinae and grouped into 
the Gamma-CoV genus (Carstens, 2010). The AvCoV genome con¬ 
sists of positive-sense single-stranded RNA with an approximate 
size of 27 kb (Jackwood et al., 2012). The genome is subdivided 
into nine open reading frames (ORFs) that encode a set of non- 
structural proteins known as the transcriptase-replicase complex, 
which functions as the RNA polymerase (Pol-la, -lb) and four main 
structural proteins (spike surface glycoprotein (S), membrane (M), 
envelope (E), and nucleocapsid (N)) (Fauquet et al., 2005; Mardani 
et al., 2008; Mo et al., 2013). 
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Among the AvCoVs, we highlight Avian infectious bronchitis 
virus (AIBV), the aetiological agent of infectious bronchitis, which 
is one of the most important viral diseases in countries with inten¬ 
sive poultry industries (Cavanagh and Gelb, 2008; Jones, 2010), 
as well as AvCoVs (IBV-like) that were recently detected in sev¬ 
eral bird orders, such as Galliformes, Anseriformes, Columbiformes, 
Charadriiformes, Passeriformes, Pelacaniformes, Ciconiiformes and 
Psittaciformes, which may influence strongly AvCoV epidemiology 
(Domanska-Blicharz et al., 2014). 

Due to its highly mutagenic character, the control of AvCoVs 
is a constant challenge for the poultry industry. The evolutionary 
dynamics of AvCoVs are characterised by the continuous emer¬ 
gence of new variants, resulting in multiple serotypes that do not 
confer cross-protective immunity to each other (Jackwood et al., 
2012). AvCoV diversity is mainly driven by three factors: the ele¬ 
vated mutation rate of the RNA-dependent RNA polymerase; the 
great diversity of birds, which provides many potential target hosts 
for different AvCoV genotypes; and the facilitation of viral transmis¬ 
sion and dissemination across large geographic areas and between 
different species by bird migration and flock lifestyles (Keesing 
et al., 2010). These factors have driven the generation of a wide 
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variety of AvCoV variants, resulting in persistent reservoirs that 
can become relevant for human populations when wild animals 
acquire synanthropic habits and start to live in close proxim¬ 
ity to domesticated poultry (Jourdain et al., 2007; Felippe et al., 
2010 ). 

Despite the health and economic importance of AvCoVs to the 
poultry industry, knowledge about the prevalence and distribution 
of these viruses in wild and synanthropic birds remains limited. 
The variety of potentially infected bird species and the annual use 
of intercontinental territory by migratory birds increase the risk of 
the introduction and spread of new pathogens and new variants 
of CoV (Jourdain et al., 2007; Somveille et al., 2013), requiring the 
development of new analytical methods. 

Phylogenetic inferences are extremely useful tools for the study 
of the molecular ecoepidemiology of pathogens and for the eluci¬ 
dation of epidemiologically important characteristics. Phylogenetic 
inferences enable the molecular characterisation of existing vari¬ 
ants and the discovery of new viral variants with potential health 
repercussions. Furthermore, these techniques permit an analysis of 
the evolutionary history of new variants (Lemey et al., 2009; Yang 
and Rannala, 2012). 

In addition to providing information on the timeline of evolu¬ 
tionary processes, phylogenetic studies can be used to compare 
and associate individuals from specific geographical regions with 
the time of the most recent common ancestor (tMRCA) of each 
region studied. Tools such as phylogeography provide a useful 
approach for understanding the origins and distributions of dif¬ 
ferent viral strains to elucidate population evolution dynamics and 
their dependence on the environment. Phylogeographic techniques 
have been applied to many viral infections that threaten human 
health, including dengue, rabies, influenza and HIV (Holmes, 2004; 
Kuhnert et al., 2011), but there is a lack of phylogeographic infor¬ 
mation regarding CoVs. 

This pilot study applied a phylogeographic approach to the 
evolution of the AvCoV spike protein gene to study the spatial 
arrangements of different genetic lineages within species and 
between closely related species, to analyse the principles and pro¬ 
cesses that determine the geographic distribution of lineages, and 
to establish a resource for improved molecular characterisation of 
available and as-yet unidentified AvCoV strains. 


2. Material and methods 

2.1. Ethical statement 

The capture and management of wild and synanthropic ani¬ 
mals for the collection of biological samples was performed in 
accordance with current Authorization and Information System of 
Biodiversity (SISBIO) recommendations, Brazil, protocol no. 34751- 
1 / 12 . 


2.2. Samples and population study 

A total of 500 tracheal and/or cloacal swabs were collected from 
42 different bird species belonging to 13 orders (Supplementary 
Table SI). The biological samples were collected in the cities of 
Campinas and Sao Paulo in the State of Sao Paulo (Southeastern 
Brazil) and the city of Florianopolis (Southern Brazil) in the State of 
Santa Catarina between 2006 and 2013. The swabs were placed in 
1,5-mL microtubes and identified with the name of the bird species 
and a specific code. The samples were immediately transported in 
coolers with ice to the Laboratory of Animal Virology (LVA) at the 
Institute of Biology at UN1CAMP and stored at -80 °C until further 
processing. 


2.3. Route map of migratory and non-migratory birds 

The maps showing the birds and their specific routes were 
derived from the BirdLife International (Cambridge, United King¬ 
dom) and NatureServe (Arlington, United States of America) 
databases. These maps consist of Geographic Information System 
(CIS) shapefiles that show the geographic extent and distribution 
of all bird species investigated in the study (Fig. 1). 

2.4. RNA extraction and detection of AvCoVs by qRT-PCR 

A total of 500 p,L of Minimum Essential Medium (MEM) (Nutri- 
cell, Brazil) was added to the microtube containing the swab with 
the biological sample. The remaining procedure was performed 
according to the QIAamp® Viral RNA Mini Kit (QIAGEN, CA, USA) 
manufacturer’s instructions. 

Each extracted RNA sample was subjected to Real Time RT-qPCR 
using previously published primers for the 5'-Untranslated Region 
(UTR) of AvCoV (Callison et al., 2006) with the addition of a double 
Quencher. Amplification of the 143-bp fragment was performed in 
a 7500 Real Time PCR Cycler (Applied Biosystems™, CA, USA) using 
the following primers and probes (IDT, Iowa, USA): 5' UTR forward 
primer, 5'-GCT ITT GAG CCT AGC GTT-3'; 5' UTR reverse primer, 
5'-GCC ATG TTG TCA CTG TCT ATT-3'; 5' UTR probe, 5'-/5HEX/CA 
CCA CCA G/ZEN/A ACC TGT CAC CTC/3IABkFQ/-3'. 

2.5. cDNA synthesis and sequencing 

The spike glycoprotein gene region was amplified from qRT- 
PCR-positive samples. Briefly, cDNA was synthesised using a High 
Capacity cDNA Reverse Transcription kit (Applied Biosystems™, 
CA, USA). Next, a nested RT-PCR assay was performed using the 
following primers in the first step to amplify a 572-bp fragment: 
forward primer S7 5'-TACTACTACCAGAGTGC(C/T)TT-3' and reverse 
primer S6 5'-ACATC(T/A)TGTGCGGTGCCATT-3'. The following 
primers were used in the second step to amplify a 530-bp fragment: 
forward primer S9 (5'-ATGGTTGGCATTT(A/G)CA(C/T)GG-3') and 
reverse primer S5 (5'-GTGCCATTGACAAAATAAGC-3') (Bochkov 
et al., 2007). This assay was performed in an Eppendorf 
MasterCycler®Pro S (NY, USA). 

The amplified PCR products were purified using the MiniElute 
PCR Purification kit (QIAGEN, CA, USA) and quantified using an 
Epoch spectrophotometer (Biotek, Vermont, USA). Purified PCR 
products were sequenced by the Central Laboratory of High Perfor¬ 
mance Technologies in Life Sciences (LaCTAD/UNICAMP, Sao Paulo, 
Brazil). Sequencing reactions were performed using the Big Dye® 
Cycle Sequencing Terminator 3.1 kit (Applied Biosystems™, CA, 
USA) following the instructions provided by the manufacturer. The 
3730XL DNA Analyser (Applied Biosystems™, CA, USA) platform 
was used for sequencing. 

2.6. Sequence data collection and phylogenetic analyses 

All bioinformatics and phylogeny experiments were performed 
using the Linux Cluster from the High Performance Computer 
Center, University of Florida, USA (UF HPC). For the bioinformat¬ 
ics and phylogenetic analyses, all available AvCoV spike 1 gene 
(SI) sequences from different animals were downloaded from the 
GenBank database (http://www.ncbi.nlm.nih.gov/genbank/) (Sup¬ 
plementary Table S2) and filtered using the following criteria: 
known sampling year; unknown sampling year; and known samp¬ 
ling year and isolation source (Fig. 2). Spike 2 gene regions (S2), 
clones, patents, and any recombinant sequences were excluded. 
Additionally, different CoV sequences were downloaded to rep¬ 
resent CoVs from the genera Alpha, Beta, Gamma, and Delta. The 
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Fig. 1. Distribution map of migratory and non-migratory birds assessed in this study. The colours represent the specific region ranges in which the birds are present. The 
highlighted names indicate birds that participate in migration events. (For interpretation of the references to colour in this figure legend, the reader is referred to the web 
version of this article). 


alignment was conducted using the parallel implementation of 
software ClustalW-MPI (Message Passing Interface) (Li, 2003). 

The initial phylogenetic analyses were performed using Max¬ 
imum Likelihood (ML) implemented in RAxML v.8 (Stamatakis, 
2014) and FastTree v.2.1.7 (Price et al., 2010) software and the 
Relaxed Neighbour-Joining (RNJ) method in Clearcut software 
using the Kimura correction for pairwise distances (Sheneman 
et al., 2006). The ML analysis in RAxML was conducted using 
the general time reversible (GTR) + G and in FastTree using the 
standard implementation GTR+CAT with 20-parameter gamma 
optimisation and a mix of nearest-neighbour interchanges (NNI) 
and sub-tree-prune-regraft (SPR) for the tree topology search. ML 
tree reliability was analysed using bootstrap support (for RAxML) 
and Shimodaira-Hasegawa test (SH-like) (for FastTree) values with 
1000 replicates (Anisimova and Gascuel, 2006); both methods com¬ 
pare the likelihood values of the best and second-best alternative 
arrangements around the branch of interest (Veras et al., 2011 ). 

Two different ML tree construction software programmes 
were used to assess the potential for different clustering pat¬ 
terns of the sequences sampled from this study with the AvCoV 
SI sequences downloaded from the GenBank database. Figtree 
v.1.4.2 (http://tree.bio.ed.ac.uk/software/figtree/) and OneZoom 
(Rosindell and Harmon, 2012) were used to visualise the large 
ML tree (file in the Supplementary Material). After the large ML 
tree was obtained, clades with node reliability >90% were selected 
using Phylopart to select sequences with high bootstrap values for 
Bayesian analysis (Prosperi et al., 2011 ). 

2.7. Bayesian approach 

Bayesian inference was performed using BEAST v.2.1.3 software 
(Bouckaert et al., 2014) assuming a strict or relaxed (uncorrelated 
Exponential and Log-Normal) molecular clock with different coa¬ 
lesced tree priors (constant population size, exponential growth, 
and Bayesian Skyline Plot (BSP)) (Drummond et al., 2005). Markov 


Chain Monte Carlo (MCMC) was executed with 450 million gen¬ 
erations and sampled every 45,000 MCMC steps. Alternatively, 
LogCombiner v.2.1.3 (available in the BEAST v.2.1.3 package; 
http://beast.bio.ed.ac.uk/logcombiner) was used to resample the 
posterior distribution states at a lower frequency. The BEAST output 
files were analysed in Tracer v.l .6 with 10% burn-in (Rambaut et al., 

2013) . The Effective Sample Size (ESS) was >500 for each param¬ 
eter assayed. Next, TreeAnnotator v.2.1.2 software (available in 
the BEAST v.2.1.3 package; http://beast.bio.ed.ac.uk/treeannotator) 
was used to produce maximum clade credibility (MCC) trees from 
each coalesced parameter. The trees were visualised in the graph¬ 
ical viewer Figtree v.1.4.2. Densitree software (Bouckaert et al., 

2014) was used for density and qualitative analyses of the trees. 
The evolution models implemented were indicated by jModelTest 
v.2 software (Santorum et al„ 2014). 

Reconstruct Ancestral State in Phylogenies (RASP) v.2.1 software 
was used to infer biogeographical histories through the Bayesian 
Binary MCMC method (BBM) to predict biogeographical ancestral 
ranges with multiple MCMC iterations (Yu et al., 2014). The MCMC 
chains were run for five million generations, and the state was sam¬ 
pled every 100 generations. The analyses were performed under the 
Fixed JC + G (Jukes-Cantor + Gamma) model and the BBM with null 
root distributions. 

2.8. Molecular clock analysis 

The tMRCA and the evolutionary rate (assuming the nucleotide 
substitution per site per year) were inferred in the sequences des¬ 
ignated as Clade 1 (representing wild and synanthropic birds) and 
Clade 2 (representing chickens) using BEAST v.2.1.3. Each sequence 
from Clades 1 and 2 was selected for the runs assuming the strict 
and relaxed molecular clock as described above. The nucleotide 
substitution models implemented in jModelTest v.2 were the 
GTR (Clade 1) and HICY (Clade 2). The relaxed (for Clade 1) and 
strict (for Clade 2) molecular clocks with the coalescent constant, 
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Fig. 2. AvCoV SI gene sequences retrieved from GenBank and classified by (a) known sampling year, (b) unknown sampling year, and (c) known sampling year and isolation 
source. 


exponential populations and BSP priors were used for Clade 1 
analysis, while the BSP coalescent prior was used for Clade 2 analy¬ 
sis. The MCMC was performed with 450 million generations and 
sampled every 45,000 MCMC steps. The ESS was >500 for each 
parameter assayed to permit sufficient mixing between the Markov 
chain. Tracer vl.6 software was used to visualise the results. 

2.9. Parsimony analysis 

To trace the history of reconstruction of character states at 
ancestral nodes, a parsimony analysis was performed using the 
Mesquite v.2.75 software package (Maddison and Maddison, 2011). 
These analyses enabled the interpretation of the evolution of the 
various characteristics of the AvCoV strains and their evolution¬ 
ary behaviour given a group of AvCoV organisms (Supplementary 
Fig. SI). 

2. 10. Bayesian phylogeography analysis 

Evolutionary history parameters, including the evolution time of 
the AvCoV strains and their spatial and temporal dispersion, were 
estimated using the Bayesian phylogeography framework imple¬ 
mented in the BEAST v.2.1.3 software package (Bouckaert et al., 


2014). To adjust the time scale of the trees, sequence data were 
obtained from the GenBank database and organised by specific 
criteria (see item 2.6). To accommodate variation between strains, 
a molecular clock approach (described above) was used to esti¬ 
mate ancestral spatial locations along the phylogenetic history. The 
results were analysed using the Tracer v.1.6 and TreeAnnotator 
v.2.1.2 software. 


3. Results and discussion 

3.1. AvCoV prevalence and phylogenetic analysis 

Samples collected from 42 species of birds were tested for the 
presence of the 5' UTR of AvCoV by real time RT-qPCR. Of these, 
65 samples were positive (13.0%), including 43 (8.6%) tracheal and 
22 (4.4%) cloacal swabs. AvCoV was detected in 22 different bird 
species distributed in 10 orders and 13 bird families. Among bird 
families, the virus was detected in Anatidae (3 of 10 species), Ardei- 
dae (1 of 3 species), Cathartidae (1 of 1 species), Columbidae (3 of 
3 species), Falconidae (1 of 2 species), Icteridae (1 of 1 species), 
Passeridae (1 of 1 species), Phalacrocoracidae (1 of 1 species), Pici- 
dae (1 of 1 species), Psittacidae (3 of 5 species), Ramphastidae (1 of 
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Fig. 3. ML phylogenetic trees of AvCoV SI gene sequences from AvCoVs detected in different animals. Phylogenetic inference using 2079 AvCoV SI gene sequences; each 
colour represents different species of birds, and the arrows indicate the positions of the birds in the tree (a). Phylogenetic inference from the large ML tree highlighting the 
sequences obtained in the study in greater detail (b). The Shimodaira-Hasegawa (SH-like) branch supports are shown on the trees. Only values greater than or equal to 75% 
are shown. Abbreviations: Pon, Pigeon; Vul, Vulture; Spa, Sparrow-Hawk; Owb, Owl bush; Owl, Owl; Wpe, Woodpecker; Gre, Great kiskadee; Par, Parrot; Red, Red-ruffed 
fruitcrow; Haw, Hawk; BR, Brazil. Only values greater than or equal to 75% are shown. (For interpretation of the references to colour in this figure legend, the reader is referred 
to the web version of this article). 


1 species), Strigidae (3 of 3 species), and Tyrannidae (2 of 2 species) 
(Supplementary Table S3). 

These findings reveal a high prevalence of AvCoVs circulating in 
birds other than chickens and suggest a potential role for wild and 
synanthropic birds as disseminators of AvCoVs in different avifau¬ 
nas. These findings demonstrate that AvCoV is widespread in birds 
from the regions studied. Considering the economic and health 
impacts of AvCoV infections, the study of these viral agents has 
important implications for epidemiology as well as the elucidation 
of molecular ecology. 


After detection of AvCoV and sequencing analysis, molecular 
evolution methods were used to phylogenetically characterise all 
positive samples. Analysis using an ML tree generated from 2079 SI 
AvCoV sequences collected from pigeon and duck, guinea fowl and 
turkey sampled between 1956 and 2014 revealed a close phylo¬ 
genetic relationship. By contrast, AvCoV sequences collected from 
wild birds in this study exhibited great divergence when com¬ 
pared with those detected in turkeys (Fig. 3). The same phylogenetic 
approach using the relaxed neighbour-joining method and an ML 
tree generated using RAxML software revealed similar patterns 
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Fig. 4. Molecular phylogenetic analysis by the ML method with AvCoV SI gene sequences from the Alpha, Beta, Gamma, and Delta CoV genera. Evolutionary history was 
inferred using the ML method based on the Kimura 2-parameter model. The black cycles represent the sequences collected in the study. The numbers above the branches 
indicate the bootstrap branch supports of 1000 replicates. Only values greater than or equal to 75% are shown. Abbreviations: Pon, Pigeon; Vul, Vulture; Spa, Sparrow-Hawk; 
Owb, Owl bush; Owl, Owl; Wpe, Woodpecker; Gre, Great kiskadee; Par, Parrot; Red, Red-ruffed fruitcrow; Haw, Hawk; BR, Brazil. 
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Fig. 5. ML phylogenetic tree of sequences from two well-supported clades extracted from the tree analysed in Fig. 3. Clade 1 (representing wild and synanthropic birds) and 
Clade 2 (representing chickens) are highlighted. The legends show the abbreviations of the animals and countries evaluated in the study. The numbers in the tree represent 
the Shimodaira-Hasegawa (SH-like) branch supports. Only values greater than or equal to 75% are shown. 


of divergence and dose relationships between species (data not 
shown). Then, a phylogenetic tree with the spike sequences from 
the Alpha, Beta, Gamma and Delta CoV groups was constructed 
to evaluate the representativity of the sequences from the study 
compared to the CoV sequences belonging to other CoV groups 
(Fig. 4). 


Two well-supported clades was detected in the large ML tree. 
The first clade was designated Clade 1 and represents wild and 
synanthropic bird sequences (including two chicken sequences: 
one from Sweden (1999) (accession number JN022543.1) and one 
from Thailand (2008) (accession number GQ885131.1)). The sec¬ 
ond was designated Clade 2 and represents chicken sequences, 
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Fig. 6. MCC trees from Bayesian framework analysis of the AvCoV SI gene. Data for the relaxed molecular clock for (a) coalescent constant population and (b) coalescent 
exponential population priors are shown. AvCoV sequences were collected from different animals. Obs.: Ckn (Sweden) indicates A1BV detected in a chicken from Sweden. 
The years preceding the branches indicate the collection date of the sequences. The posterior probability is represented by the numbers in the branches. 


with the exception of a duck sequence from China collected in 2004 
(accession number JF699751.1) (Fig. 5 and Supplementary Fig. S2). 

3.2. Tree topologies based on density and inference of 
biogeographic history (RASP) 

Because two clades were detected, a Bayesian approach using 
Densitree software was applied to verify the tree topologies based 
on the density of the stored trees (10,000 trees) from BEAST anal¬ 
ysis. In this analysis, the branch with the strongest and densest 
colour was assumed to be the most likely tree topology. A chicken 
sequence from Sweden (1999) was observed as an outgroup in both 
coalescent priors models used (See Supplementary Fig. S3). 

Based on the results of the Densitree software analysis described 
above, MCC trees were constructed by imposing a relaxed molec¬ 
ular clock. An A1BV strain from Sweden (1999) was inferred as the 
most recent common ancestor of the sequences collected in this 


study (Fig. 6 and Supplementary Fig. S4), suggesting that a strain 
collected 15 year ago may have been responsible for the emergence 
of new AvCoV strains in wild and synanthropic animals. 

Following these analyses, the BBM framework was imple¬ 
mented to obtain a biogeographical reconstruction of the ancestral 
ranges of the AvCoVs detected in this study. First, ancestral nodes 
of A1BV from chickens were analysed. Ancestral nodes 78 (C, D, E, 
and F indicate a European ancestor), 80 (I, J, and L have a ancestral 
range from Australia, Canada, and South Korea), 92 (P and N have 
an ancestral range from Iraq and China) and 93 (O, P, and Q have 
an ancestral range from Slovenia, Italy, and Iraq) exhibit ambiguity 
and suggest the presence of two, three and four different recon¬ 
struction ancestral ranges (i.e., for each node analysed, there exists 
more than one biogeographic ancestor). In other words, the BBM 
analyses inferred two, three and four dispersal events involving 
countries on different continents (North America, Europe, Oceania, 
and Asia) (Fig. 7). 
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Fig. 7. Biogeographical ancestor reconstructions implemented in RASP software using the BBM method. Graphical results for the ancestral distributions at each node obtained 
from the phylogenetic analyses of AvCoV strains by BBM analysis, (a) RASP area reconstructions, (b) Each colour followed by a letter abbreviation represents a specific country 
assuming the possible ancestral ranges for the different nodes. The blue traces in the trees are the sequences collected in this study. The probabilities of alternative ancestral 
ranges are shown in pie charts at each node, (c) Possible route of the AvCoV strain from the host bird to the Brazilian birds assessed in the study, (d) Results from ancestral 
distributions at node 111 showing the percentage of ancestral identity between Swedish and Brazilian virus. Yellow arrows highlight the results from 111 and 109 ancestral 
nodes. Abbreviations: BR, Brazil; SE, Sweden; TH, Thailand. (For interpretation of the references to colour in this figure legend, the reader is referred to the web version of 
this article). 
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Fig. 8. BSP reconstruction of AvCoV strains detected in chickens (a) and in wild and synanthropic birds (b). Purple branches, chicken sequences; Yellow branches, pigeon 
sequences; Blue branches, sequences from wild birds. N, number of transmissions between different bird species. The blue arrow represents the increase in transmission 
between different bird species. The effective number of the population is represented on the Y-axis. Time is represented on the X-axis. Black lines indicate the medians and 
blue lines the credibility based on 95% highest posterior density (HPD) intervals. (For interpretation of the references to colour in this figure legend, the reader is referred to 
the web version of this article). 
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BBM analyses of the ancestral ranges of AvCoV strains from wild 
and synanthropic birds (highlighted with blue traces) revealed that 
most strains had only Brazilian identity. However, ancestral nodes 
109 (A and R) and 111 (A and D) had a biogeographic ancestrality 
range from other countries. Data from the MCC trees showed that 
an AIBV strain from Sweden detected in chickens was the nearest 
ancestor; based on Node 111, the ancestrality frequency was 0.84% 
(AD), 0.19% (D), and 0.64% (from another ancestor range) (Fig. 7). 
From these analyses, our results suggest that wild birds may play a 
role in the spreading of AvCoVs amongst poultry around the globe. 

3.3. Molecular clock analysis 

The molecular clock hypothesis was evaluated for AvCoV strains 
using the likelihood ratio test (LRT) (Posada and Crandall, 1998), 
which compares the likelihood values of the topology assuming 
and not assuming the molecular clock. The time of divergence of 
the sequences from Clade 1 (AvCoV detected in wild and synan¬ 
thropic birds) and 2 (AIBV sampled from chickens) was estimated, 
and the evolutionary rates for the sequences from Clade 1 and Clade 
2 were calculated at 1.59 x 10 -2 and 2.21 x 10~ 4 substitutions per 
site per year, respectively (Supplementary Table S4). The sequences 
from Clade 2 followed the molecular clock model according to 
the CoV evolutionary rate (Jackwood et al., 2012). Previous stud¬ 
ies have estimated an evolutionary rate of 0.3-0.6 x 10~ 2 for some 
AIBV strains (Cavanagh et al., 1998). Because the sequences were 
obtained from wild and synanthropic birds and AIBV is primar¬ 
ily found in chickens, we would expect to obtain a different result 
(i.e., strains from different environments may exhibit differences 
in evolutionary rates). 

A high rate of CoV evolution leads not only to the generation of 
a diversity of strains and genotypes of CoV species but also to new 
species that can adapt to new hosts and ecological niches (Woo 
et al., 2009). Recent studies have identified CoV isolates from pea¬ 
cock ( Pavo cristatus ), partridge ( Alectoris sp.), Guinea fowl ( Numida 
meleagris ) and teal ( Anas sp.) with AIBV-like genomic characteris¬ 
tics (Liu et al., 2005; Chu et al., 2011). In addition, close contact 
between infected and susceptible animals facilitates the transmis¬ 
sion of RNA viruses between species when two hosts are closely 
evolutionarily related (Cavanagh, 2005), as is the case for chickens 
and quails, which belong to the Phasianidae family. 

3.4. Reconstruction of the demographic history 

To reconstruct the demographic history and analyse past popu¬ 
lation growth over time, the BSP method was applied by enforcing 
the relaxed (Clade 1 sequences, AvCoV detected in wild and synan¬ 
thropic birds) and strict (Clade 2 sequences, AIBV sampled from 
chickens) molecular clocks. This analysis demonstrated a constant 
and stationary equilibrium in demographic chicken populations 
over time (1956-2014, sequences from Clade 2) (Fig. 8A). Inter¬ 
estingly, the BSP reconstruction of the demographic history of 
AvCoV revealed an exponential increase in effective population 
size in 2009 (sequences from Clade 1) (Fig. 8B). An MCC tree from 
the BSP prior was positioned with the BSP graphic to verify data 
reproducibility. The results suggested that the potential increase of 
AvCoV strains in wild birds was directly related to the increase in 
animal interspecies AvCoV transmission (n = l to n = 12 to n = 32) 
(Fig. 8B). Use of the coalescent parameter from the exponential 
population growth model to observe the direction of viral trans- 
missibility using the MCC tree approach (chicken (n = 1) to pigeon 
(n = 12) to wild birds (n = 32)) enabled an estimation of a coherent 
epidemiological trace for interspecies AvCoV transmission (Fig. SB). 

The frequency of viral transmission from birds to different hosts 
is increasing, and the viruses are gaining the ability to infect ani¬ 
mals of different species, including the human population (e.g., 


the H1N1, H2N2, and H3N2 pandemic influenza viruses), either 
by viral evolutionary mechanisms that result in a new viral vari¬ 
ant or via contact with different hosts. The main sources of new 
human viral diseases are enzootic and epizootic animal viruses 
(Parrish et al., 2008), and thus an in-depth phylogenetic and phy- 
logeographic examination of the epidemiological context of AvCoV 
genetic diversity is extremely important. 

This study employed an in-depth phylogenetic approach, and 
the results may be extended not only to understand the molecu¬ 
lar dynamics of AvCoV evolution in different animals but also to 
other epidemiologically important pathogens. The emergence of 
new viral variants requires a better understanding of the molecular 
eco-epidemiology of these viral agents, which will require addi¬ 
tional data regarding viral evolutionary behaviour, given the broad 
transmission and dissemination ranges of AvCoVs and CoVs in gen¬ 
eral. Finally, due the shortage of scientific studies applying in-depth 
bioinformatics and phylogenetic approaches in the AvCoV field and 
considering that many studies that detected CoVs in wild birds 
were conducted in recent years, a sudden expansion of AvCoVs has 
emerged in the public domain. Thus, some of the conclusions might 
be due to sampling bias. 
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